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We  present  new  convergence  and  stability  results  for  least  squares  inverse  problems  in¬ 
volving  systems  described  by  analytic  semigroups.  The  practical  importance  of  these  results 
is  demonstrated  by  application  to  several  examples  from  problems  of  estimation  of  material 
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1^  Introduction 

In  this  presentation  we  consider  abstract  inverse  problems  in  a  least  squares  formulation 
for  parameter  dependent  partial  differential  equations.  We  are  interested  in  approximation 
ideas  which  lead  to  viable  computational  techniques  for  such  problems.  We  pursue  our 
investigations  in  the  context  of  the  general  framework  for  convergence  and  stability  developed 
by  Banks  and  Ito  im[BI).  Motivated  by  questions  related  to  the  use  of  accelerometer  data  to 
estimate  parameters  in  flexible  structures,  we  focus  on  second  order  (in  time)  systems  with 
sufficient  damping  so  that  the  system  can  be  modeled  by  an  analytic  semigroup. 

We  state  and  prove  a  new  approximation  result  (a  Trotter-Kato  type  theorem)  for  ana¬ 
lytic  semigroups.  This  theorem  gives  conditions  under  which  a  family  of  approximating  semi¬ 
groups  and  all  its  time  derivatives  converges  to  a  limit  semigroup  and  all  its  time  derivatives, 
respectively.  These  theoretical  results  are  then  stated  in  terms  of  simple,  readily  checked  con¬ 
ditions  on  the  sesquilinear  forms  defining  “stiffness”  and  “damping”  in  the  abstract  second 
order  systems. 

We  discuss  several  examples  which  indicate  clearly  the  practical  importance  of  these  new 
convergence  results.  In  two  of  the  examples  presented  (the  damped,  cantilevered  Euler- 
Bernoulli  beam  and  the  “RPL  experiment”  structure),  we  can  apply  the  abstract  theorem 
to  substantially  sharpen  results  already  found  in  the  research  literature.  In  a  third  example 
involving  a  two  dimensional  grid  structure,  we  note  that  the  theory  can  be  applied  to  obtain 
new  results  for  acceleration  data  convergence  in  the  least  squares  inverse  problems. 

2.  Abstract  Inverse  Problems 

In  this  section  we  formulate  a  class  of  inverse  problems  as  abstract  least  squares  optimiza¬ 
tion  problems  constrained  by  evolution  equations  in  a  Hilbert  space  and  summarize  some  of 
our  previous  results  for  such  problems  [B],  [BI],  [BK]. 

We  assume  we  are  given  a  parameter  dependent  system 

(9  1  x  u(t )  =  A(q)u{t)  +  F{t,  q),  0  <t<T, 

1  ]  «(0)  =  u0(q) 

for  states  tt(i;  q)  in  a  Hilbert  space  H.  The  parameters  q  are  to  be  chosen  from  an  admissible 
parameter  set  Q  contained  in  a  metric  space  (Qi,d).  We  assume  throughout  that  Q  is  a 
compact  subset  of  Q\  (a  “regularization”  assumption). 

We  are  given  observations  {zj},  a  set  of  points  in  the  observation  space  Z,  along  with  an 
observation  map  C  :  C(0,T]H)  — >  Z  from  the  states  to  the  observations.  The  points  Z{  are 
observations  for  Cu(ti\  q),  0  <U  <  T,  and  it  is  this  “data”  to  which  we  wish  to  fit  the  model 
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by  a  best  choice  of  the  parameter  q.  Formally,  the  problem  can  be  stated  as 


(?) 


Find  q  €  Q  to  minimize  over  Q  the 
functional  J{q)  =  Yi  |Cu(*»!  <?)  —  Zi\z- 


The  observation  operator  C  is  of  fundamental  importance  to  the  discussions  in  this  pa¬ 
per.  For  parabolic  systems  (2.1),  a  typical  example  for  C  arises  from  pointwise  evaluation  in 
the  spatial  variables,  e.g.,  Cu{t\q)  —  {u(t,  Xj]  g)}*=1  where  Z  =  Rl.  For  problems  involving 
structures  such  as  beams  or  plates,  several  examples  arise  in  practice.  If  one  takes  measure¬ 
ments  with  a  laser  vibrometer,  then  one  obtains  measurements  at  specific  points  in  space  for 
the  velocities  ut(U ,  Vj\  q)  so  that  the  map  C  is  a  composite  of  time  differentiation  followed 
by  pointwise  evaluation.  If  the  measuring  devices  are  accelerometers,  one  has  observations 
for  *be  accelerations  utt(Li,Xj4,q)  and  thus  C  is  time  differentiation  (twice)  followed  by 
pointwise  evaluation. 

The  inverse  problems  outlined  here  are  generally  infinite  dimensional  in  both  the  states 
u(t)  €  H  and  the  parameters  q  €  Q  and,  moreover,  involve  unbounded  operators  related  to 
the  states  (-A(g))  and  the  observations  ( C ).  Thus  to  develop  efficient  computational  methods, 
one  must  make  finite  dimensional  approximations  for  both  the  state  and  parameter  spaces, 
H  and  Q}  respectively.  For  the  discussions  of  this  paper,  we  shall  restrict  our  considera¬ 
tions  to  approximations  of  the  state  space.  In  some  cases  (where  the  parameter  sets  Q  are 
finite  dimensional,  either  naturally  or  through  some  a  priori  parameterization)  no  approxi¬ 
mation  of  the  parameter  set  is  necessary.  Methods  for  approximation  of  Q  by  sets  QM  have 
been  fully  discussed  elsewhere  (see  [BK])  and  these  ideas  could  readily  be  incorporated  in 
our  presentation.  Since  this  would  add  nothing  to  either  the  difficulties  or  their  resolution 
that  are  the  focus  here,  to  avoid  considerable  notational  tedium  we  do  not  consider  such 
approximation  ideas  in  the  discussions  below. 

Thus,  we  consider  approximate  (to  be  made  precise  below)  state  spaces  HN  C  H  with 
associated  approximate  states  u N{t',q)  €  HN  satisfying  the  approximate  systems 

,  *"(<)  =  A"(,)u’‘(t)  +  PNF(t,<l) 

(2'2)  u"(0)  =  0(9) 


where  AN  is  an  approximation  to  A  and  PN  is  the  orthogonal  projection  of  H  onto  HN.  The 
corresponding  approximate  inverse  problems  are  given  by 


( VN ) 


Find  (f*  €  Q  to  minimize  over  Q  the 
functional  JN(q )  =  Yj  |Cu^(t,;  q )  -  zi\z- 


Of  course,  there  are  a  number  of  reasonable  ways  in  which  the  approximations  may  be 
made  to  arrive  at  the  problems  ( VN ).  In  analyzing  different  methods  and  their  behavior 
in  the  context  of  inverse  problems,  a  number  of  questions  related  to  parameter  convergence 
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arise  naturally.  For  example,  given  a  fixed  set  of  data  z  =  {z*},  do  optimal  parameters 
(f*  of  (VN)  converge  in  some  sense  to  an  optimal  parameter  q  for  (' V )?  More  generally,  one 
might  also  incorporate  the  continuous  dependence  of  the  estimates  on  the  observations  in  the 
concept  of  method  stability  [B],  [BK].  For  this,  one  denotes  by  qN(zK)  optimal  parameters 
obtained  from  (' PN )  for  observations  zK  and  requires  that  qN{zK)  — >  g(z°)  in  some  sense  as 
N  — *  oo  and  zK  — >  z°  in  Z,  where  g(z°)  is  an  optimal  parameter  for  (V)  corresponding  to 
observations  z°.  These  issues  are  carefully  discussed  in  [B],  [BK]  where  it  is  shown  that  for 
both  parameter  convergence  and  method  stability,  it  suffices  to  argue  that 

For  arbitrary  sequences  {g^}  in  Q  converging 
(2.3)  to  g  £  Q,  we  have  CuN(t\ qN)  — >  Cu{t\  q)  as  N  — *  oo 

for  each  t  £  (0,  T ). 

Thus,  certain  fundamental  aspects  of  approximation  in  inverse  problems  can  be  reduced  to 
the  convergence  statement  in  (2.3)  and  we  shall  deal  with  conditions  under  which  (2.3)  can 
be  guaranteed  in  the  subsequent  discussions  of  this  paper. 

In  [BI],  Banks  and  Ito  developed  a  general  framework  for  convergence  and  stability  which 
we  summarize  here  in  an  extended  form  to  allow  treatment  of  complex  valued  operators  in 
complex  state  spaces  (the  sesquilinear  forms  introduced  in  [BI]  were  tacitly  assumed  to  be 
generated  by  partial  differential  equations  with  real  coefficients  and  the  theory  developed 
there  was  adequate  for  such  examples).  As  in  [BI],  our  goal  is  to  isolate  properties  of  the 
system  (2.1)  and  the  approximation  families  to  insure  that  (2.3)  holds.  We  state  these 
properties  in  terms  of  conditions  on  a  parameter  dependent  sesquilinear  form  arising  from 
reformulating  (2.1)  in  a  weak  or  variational  sense.  To  this  end  we  consider  the  system  as  a 
variational  equation  holding  for  all  cf>  £  V  : 

r9  4\  (u(t),<£)  +  or(g)(u(t),<£)  =<  F(t,  q),<f>  >, 

K  *  *(0)  =  u0(q). 

Here  V  and  H  are  complex  Hilbert  spaces  with  V  imbedded  densely  and  continuously  in 
the  pivot  space  H\  i.e.  V  <— *  H  =  H*  <— >  V*.  The  sesquilinear  form  a(q)  :  V  x  V  — *  C  is 
related  to  the  operator  A(g)  as  described  below,  <  .  >  is  the  inner  product  in  H,  (•,•)  is  the 
duality  product  (-jOv^.v  and  equation  (2.4)  is  interpreted  in  the  V*  sense.  The  conditions 
on  a  needed  are: 

(A)  Continuity  in  g:  For  g,  g  £  Q,  we  have  for  all  <p,  ip  £  V 

K?)(^V0  -  r(g)(^,VOI  ^  d(9»?)l 4>\vH\v- 

(B)  V-Coercivity:  There  exists  c\  >  0  and  some  real  A0  such  that  for  g  £  Q,  $  £  V,  we  have 

Re  &(q)(<t>,  <f>)  +  AqI^Ih  >  Cj \<f>\^. 
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(C)  Boundedness  in  q:  There  exists  c2  >  0  such  that,  for  q  E  Q,  E  V,  we  have 

W(q)(<f>,ip)\  <  c2|^jv|^!v- 

Several  remarks  concerning  these  conditions  are  useful.  First,  we  can  weaken  (A)  slightly 
(see  [BRR]):  From  (C)  it  follows  that  one  can  define  A(q )  in  £(V,V*)  by  a(q)(<t>,ip)  = 
and  then  one  can  replace  condition  (A)  by  the  requirement  that  q  — *  A(q) 
be  continuous  from  Q  to  C(V,  V*). 

Moreover,  for  a  :  V  x  V  — *  C  defined  using  real  coefficients,  i.e.  real  scalars,  it  suffices 
to  replace  (B)  by  the  condition:  <r(q)(<i>,<f>)  +  Ao|<^|tf  >  Ci|$|y-  for  all  real  valued  functions  <f> 
in  V  (this  is  condition  (B)  as  stated  in  [BI]).  To  see  this,  we  note  that  cr(q)((j>  +  itp,  <j>  +  iip) 
=  <r(q)(4>,4>)  +  and  thns  Re  a(g)(</>  +  iip,  <f>  +  tip)  = 

a(q)(<f>,(j))  +  cr(q)(ip,i>).  Therefore,  cr(q)(<f>,<f> )  +  A0|^|^  >  Ci|^|yr  for  all  real  valued  <f>  in  V 
implies  immediately  that  condition  (B)  above  holds. 

Under  conditions  (B)  and  (C),  a  defines  operators  A(q)  such  that  cr(g)(^>,  -0)  = 
<  —A{q)<f>,ij>  >  for  <f>  E  dom (A(q)),ifr  E  V  with  dom(A(g))  dense  in  V  (e.g.,  see  [S]).  More¬ 
over,  A(q)  is  a  sectorial  operato-  with  (A I  —  A(g))dom(.A(g))  =  H  for  all  A  with  ReA  >  A0; 
indeed  R\(A(q))  =  (A I  —  A(q))~x  exists  as  a  bounded  operator  on  H  for  all  A  in  the  com¬ 
plement  of  a  sector  with  A0  as  a  vertex.  The  operator  A(q)  generates  an  analytic  semigroup 
T(t‘,q)  which  can  be  used  in  defining  mild  solutions  of  the  system  (2.4)  or  (2.1);  i.e., 

(2.5)  u(i)  q)  =  T(t ;  q)u0{q)  +  /  T(t  -  s;  q)F(s,  q)ds. 

Jo 

Galerkin  approximations  (e.g.,  systems  (2.2))  for  (2.5)  can  be  developed  in  the  context  of 
sesquilinear  forms  satisfying  conditions  (B)  and  (C).  Let  HN  be  a  family  of  finite  dimensional 
subspaces  of  H  satisfying  HN  C  V  and  the  condition 

.  .  For  each  <f>  E  V,  there  exists  <}>N  E  HN  such  that 

| <f>  —  $N\v  — >  0  as  N  — »  oo. 

(We  note  that  many  popular  approximation  schemes  -  e.g.,  linear  splines  and  more  generally, 
many  finite  element  schemes  -  satisfy  condition  (Cl).)  To  define  AN(q)  :  HN  — *  HN ,  we 
restrict  cr(q)  to  HN  x  HN  and  denote  this  restriction  by  0N{q).  Then  oN  satisfies  conditions 
(A),  (B),  (C)  on  HN  x  HN  which  implies  existence  of  a  bounded  linear  operator  AN(q) 
on  Hn  satisfying  0N{q){<f>,  t/>)  =<  —^(q)#,^  >  for  all  <f>,  rp  E  HN .  From  (B)  and  (C)  it 
follows  that  the  AN(q)  are  uniformly  dissipative  (sectorial)  and  generate  analytic  semigroups 
TN(t-q)  on  HN.  Solutions  of  the  systems  (2.2)  are  then  given  by 

T 

(2.6)  uw(iig)  =  r"(t;,)P'%,(<!)+  /  TN{t-3-,q)PNFfls,q)d3. 

Jo 
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It  can  then  be  established  that  under  conditions  (A),  (B),  (C)  and  (Cl),  we  have 
v.N(t]qN)  — »  u(t\ q)  in  V  for  arbitrary  sequences  {gN}  with  qN  — >  q  in  Q.  Observe  that 
this  immediately  yields  (2.3)  if  the  operator  C  possesses  certain  boundedness  properties. 
The  focus  of  this  note  involves  cases  where  the  operator  C  does  not  have  such  boundedness. 

The  arguments  for  the  convergence  uN(t\qN)  — ►  u(t]  q)  for  the  situation  here  (complex 
Hilbert  spaces  V,  H  and  complex  valued  sesquilinear  form  a)  are  essentially  the  same  as 
those  in  [BI],  Slight  changes  in  the  arguments  for  Theorem  2.2  of  [BI]  are  necessary  to  treat 
the  case  of  complex  valued  inner  products  and  sesquilinear  forms  (essentially  one  need  only 
use  Re  cr( •,  •)  and  Re  <  -,  •  >  in  some  of  the  inequalities).  In  these  arguments  one  does  not 
use  directly  the  analyticity  properties  of  the  semigroups  T(t;  q),TN(t]q).  Rather  one  relies 
heavily  on  resolvent  estimates  of  Tanabe  [T]  depending  on  the  V-coercivity  of  a  along  with  a 
resolvent  convergence  form  of  the  Trotter-Kato  approximation  theorem  of  linear  semigroup 
theory  (see  [BI]  for  details). 

Our  interest  here  is  mainly  in  second  order  systems  of  the  form 

u(t)  +  B(q)ii(t)  +  A(q)u(t)  =  f(t ,  q) 

(2.7)  u(0)  =  u0 

u(0)  =  v0 

where  A(q)  is  a  generalized  stiffness  operator  and  B(q)  is  a  generalized  damping  operator. 
Again  we  consider  this  equation  in  a  weak  or  variational  sense  defined  via  parameter  depen¬ 
dent  sesquilinear  forms  in  a  complex  Hilbert  space  V  ►  H  =  H*  ■— >  V* .  We  are  given  a 
stiffness  sesquilinear  form  <7i(g)  :  V  x  V  — >  C  that  is  symmetric  and  satisfies  the  boundedness 
condition  (C).  Then  there  exists  A(q)  (E  C(V,V*)  such  that  cri(q)(<f>,ip)  =  (A(q)<f>, 

We  also  have  a  damping  sesquilinear  form  ^(g)  :  V  x  V  — *  C  which  satisfies  (C)  so  that 
there  exists  B(q)  G  C(V,  V*)  with  a2 (g)(^,  i/')  =  {B{q)<f>,'^)v‘,v-  We  then  reformulate  the 
system  (2.7)  and  seek  solutions  u{t)  G  V  satisfying  for  all  $  E  V 

.  v  (u(t),  <f>)  +  cr3(q)(u(t),  <f>)  +  <7j (g)(u(t),  <f>)  =<  f{t,  q),  (f>  > 

u{0)  =  uo,u(0)  =  vo. 

As  is  standard  practice,  we  rewrite  this  in  first  order  vector  form  on  H  =  V  x  H  and 
V  =  V  x  V  in  the  coordinates  (u,  ft).  To  this  end,  define  a(q)  :  V  x  V  — ►  C  by 


(2.9)  *(«)((«, «),(& VO)  =  ~  <v,<j)>v  +<ri(q){u,'ip)  +  a2{q)(v,ip) 

so  that  (2.8)  may  be  rewritten  as 

(o  in'*  («<0.X)  +  ff(9)(w(0.x)  =<  B{t,q),x  >H 

{  U)  w{0)  =  (uo,vo) 
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for  w(t)  =  (u(t),u(t))  and  x  =  in  V  with  F(t,q)  =  (0 ,f(t,q)).  Or,  if  as  in  the  usual 

practice,  we  abuse  notation  and  do  not  distinguish  between  row  and  column  vectors,  we  may 
write  this  in  equivalent  operator  form 

w(t)  =  «4(g)u>(t)  +  F[t,q ) 

u>(0)  =  (ti0,v0) 

where  o(g)(x,£)  =  (~-4(g)x,£)  with 

0  I 
-A(q)  - B{q )  • 

For  our  treatment  here  (as  for  that  in  [BI])  we  assume  that  01(g)  satisfies  conditions 
(A),  (B),  and  (C).  We  further  assume  that  02(9)  satisfies  conditions  (A)  and  (C).  Then  the 
strength  of  the  coercivity  assumption  on  02(g)  determines  the  properties  of  the  semigroup 
generated  by  .A(g).  For  example,  if  02(g)  satisfies  (B),  then  A(q)  is  strongly  V-coercive  and 
generates  an  analytic  semigroup  on  7i.  (Actually,  .4(g)  is  V,-coercive  with  Vq  =  Vq  x  Vq 
where  is  V  taken  with  the  equivalent  inner  product  0i(g)(-,-)~  see  [BI].)  In  [BI],  it  is 
assumed  only  that  <r2  is  //-semi coercive:  There  exists  b  >  0  such  that  for  all  <f>  6  V  we  have  Re 
cr(g)(<^ ,  <f>)  >  b\(f>\2H.  In  this  case  one  can  argue  only  that  .4(g)  generates  a  strongly  continuous 
semigroup  on  H.  The  theory  for  second  order  systems  with  this  weak  damping  is  developed 
in  [BI]  in  order  to  treat  several  forms  of  damping  (spatial  hysteresis,  time  hysteresis,  bending 
rate  damping)  which  are  of  physical  interest  and  yet  do  not  satisfy  the  strong  V— coercivity 
assumption.  The  convergence  theory  obtained  yields  that  uN(t;  qN )  — *  u(t\  q)  in  V  norm, 
uN(t]  qN)  —*  ti(f;g)  in  H  norm  whenever  qN  —*  q. 

As  opposed  to  [BI],  we  wish  to  consider  in  this  paper  the  case  where  02(g)  does  satisfy 
the  strong  V— coercivity  condition  (B)  and  hence  .4(g)  is  the  infinitesimal  generator  of  an 
analytic  semigroup.  We  obtain  immediately  that  o-(g)  of  (2.9)  is  V— coercive  and  the  first 
order  theory  outlined  above  can  be  applied  directly  to  the  system  (2.10).  This  yields  the 
convergence  statement  in  V  -  i.e.,  uN(t\qN )  — ►  u[t\q)  in  V,  uN(t-,qN)  — ►  u(t\q)  in  V.  A 
natural  question  arises  as  to  whether  we  can  use  the  analyticity  of  the  semigroups  in  this  case 
to  obtain  stronger  results.  The  next  section  is  devoted  to  results  that  yield  an  affirmative 
answer  to  this  question. 

3.  Analytic  Semigroups  and  Approximation 

We  first  give  a  general  approximation  theorem  for  analytic  semigroups  that  is  a  general¬ 
ization  of  the  well-known  Trotter-Kato  theorem  [BK]. 

Theorem  3.1,  Suppose  we  have  complex  Hilbert  spaces  X  and  N  =  1,2...,  with 
XN  C  X.  Let  PN  :  X  —*  XN  denote  the  orthogonal  projection  of  X  onto  XN  satisfying 
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PN  — »  I  strongly.  Suppose  that  AN  and  A  are  the  infinitesimal  generators  of  analytic 
semigroups  SN(t)  and  S{t)  on  XN  and  X  respectively  that  satisfy  the  following: 

There  exists  a  region  53  =  53$  =  {A  £  C  :  |ar<?(A  —  A0)|  <  f  +  £}>  where  6  >  0, 

such  that  53  U{Ao}  C  p(A)  D~=1  p(AN)  and 

(i)  there  exists  a  constant  M  independent  of  N  such  that 


I*a(^)I  < 


M 

|A-Ao| 


for  all  A  €  X  an<L  N  =  1,2,...; 


(ii)  for  some  A  6  53  an d  eac h  x  £  X  we  have  Rx(AN)PNx  — >  Rx(A)x. 
Then  we  have 

(Hi)  For  each  x  £  X,  SN(t)PNx  — ►  S(t)x  uniformly  in  t  on  compact 
subintervals  o/[0, 00); 


(iv)  For  each  x  £  X  and  integer  k  >  1,  (AN)k  SN(t)PN  x  —*  AkS(t)x 
uniformly  in  t  on  compact  subintervals  o/(0,oo). 

Proof.  The  statement  of  (iii)  under  the  given  conditions  is  just  a  variant  of  the  well-known 
Trotter-Kato  theorem  and  follows  immediately  from  Theorem  II. 1.14  of  [BK].  To  argue  result 
(iv),  we  first  note  the  convergence  in  (ii)  for  some  A  £  X)  implies  that  the  convergence  holds 
for  aU  A  G  X)  •  In  light  of  the  resolvent  bounds  of  (i),  this  follows  from  the  identity  for 

/x,  A  G  X) 

Rtl{AN)PN  -  RU(A)  =  \1  +  (A  -  p.)Ru{AN)PN]\Rx{AN)PN  -  Pa(,4)][/  +  (A  -  p)R^A)} 


which  is  readily  established  using  the  standard  resolvent  identity  RX(A)  —  R^A)  =  (p  — 
X)Rx(A)Rli(A). 

The  analyticity  of  the  semigroups  SN(t )  and  S{t )  allow  us  to  write  for  t  >  0  (e.g.,  see 

IP]) 

AkS{t )  =  —  f  XkeXtRx{A)d\ 

2m  Jr 

and 

( AN)kSN(t)PN  =  f  XkeMRx{AN)PNdX. 

2m  Jr 


Here  T  is  a  positively  oriented  contour  through  A0  lying  in  53  with  arg(A  —  ^0)  =  dzv  for 
A  7^  A0  where  v  is  fixed  in  (|+6— e,  f +^).  The  desired  convergence  results  follow  immediately 
from  the  inequality 


\(AN)kSN{t)PNx  -  ,4fcS(t)x| 

~hSr  |A|‘|eA'1  "  R>(A)*\dX 
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using  the  resolvent  convergence  of  (ii)  by  noting  that  the  integrand  in  this  integral  is  domi¬ 
nated  by  an  integrable  function. 

We  now  return  to  the  second  order  systems  of  Section  2  -  see  (2.8),  (2.9),  (2.10)  -  and 
use  Theorem  3.1  to  obtain  our  main  convergence  results.  As  before  we  take  7 i  =  V  x  H  and 

V  -  V  xV.  Let  hiN  =  HN  x  HN  and  PN  be  the  orthogonal  projection  of  hi  onto  hiN . 

Theorem  3.2  Let  o\(q)  and  a2(q)  in  (2.8)  satisfy  conditions  (A),  (B)  and  (C)  and  let 
HN  C  V  satisfy  condition  (Cl).  Let  be  arbitrary  in  Q  with  qN  —>  q.  Then  we  have 

(i)  The  sesquilmear  form  o(q)  given  by  (2.9)  satisfies  conditions  (A),  ( B ),  (C)  in  the  norms  of 

V  and  hi  and  the  operator  A(q)  defined  via  a(q)(x,()  =<  —  A(q)x >n  for  x  €  dom(A(q)) 
is  the  infinitesimal  generator  of  an  analytic  semigroup  T(t\q)  on  hi. 

(ii)  Let  AN(q)  denote  the  operator  obtained  by  restricting  o(q)  to  hiN  x  hiN  and  let  TN{t\q) 
denote  the  corresponding  analytic  semigroups  on  7iN .  Then  we  have 

(a)  For  each  x  6  Ti,  TN(t\ qN)PNX  ~ ►  T(t;  q)x  in  hi  uniformly  in  t 
on  compact  subintervals  o/[0,oo); 

(b)  For  each  x  £  Tl  and  positive  integer  k,  AN (qN)kTN {t\qN)PN x  —*  A(q)kT(t;  q)x 
in  hi  uniformly  in  t  on  compact  subintervals  of  ( 0,oo). 


The  proof  of  this  theorem  is  now  rather  straightforward.  We  first  note  that  it  is  readily 
shown  with  routine  calculations  that  ai,cr2  satisfying  (A),  (B),  (C)  in  the  V  and  H  norms 
imply  that  a  satisfies  (A),  (B),  (C)  with  the  norms  of  V  =  V  x  V  and  hi  —  V  x  H.  Result 
(i)  of  the  theorem  thus  is  established. 

For  result  (ii),  we  apply  Theorem  3.1  with  X  =  hi  =  V  x  H  and  XN  =  hiN  =  HN  x  HN 
along  with  arguments  of  [BI] .  Restricting  o{qN)  to  hiN  x  hiN  to  obtain  the  operators  AN  = 
■W),  we  have  that  condition  (B)  is  satisfied  in  hiN  with  the  constants  uniform  in  N  - 
i.e.,  the  uniform  sector  condition  and  the  uniform  resolvent  bounds  in  (i)  of  Theorem  3.1  are 
readily  seen  to  hold.  For  the  resolvent  convergence  of  (ii)  we  refer  to  Theorem  2.2  of  [BI], 
noting  that  we  have  all  the  hypotheses  of  that  theorem  holding  here  in  the  sense  of  the  V 
and  hi  norms.  Hence  the  same  arguments  (modified  slightly  as  mentioned  in  Section  2  above 
to  account  for  complex  Hilbert  spaces  V  and  hi)  given  for  Theorem  2.2  of  [BI]  can  be  used 
here  to  establish  (ii)  of  Theorem  3.1.  The  convergence  statements  (a)  and  (b)  of  Theorem 
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3.2  then  follow  directly  from  (iii)  and  (iv)  of  Theorem  3.1. 

We  are  now  in  a  position  to  use  the  results  of  Theorem  3.2  to  address  the  question  of 
stronger  convergence  results  for  uN(t]  qN )  to  u(t ;  q )  raised  in  connection  with  the  observation 
operator  C  in  Section  2.  We  first  consider  the  case  of  homogeneous  systems  (i.e. ,  /  =  0 
in  (2.8))  in  which  case  wN(t]qN)  =  (uN(t-,qN),uN(t-,qN))  -  TN(t\qN)PNx  and  w(t;q)  = 
(u(t]q),u(t\q))  =  T(t\q)x ■  For  x  £  dom  (.4(g))  we  have  wN{t)  =  AN(qN)TN  (t]qN)PN  x, 
=  AN(qN)3TN(t-,qN)PNx,  etc.,  while  w(t)  =  A{q)T(t\ q)x,  w(t)  =  A{q)2T(t]  q)x, 
etc.  Thus,  from  statement  (b)  of  Theorem  3.2  we  find  ibN(t\qN)  — >  w{t\q)  in  H,  which 
yields  u^t(t\qN)  — >  utt(t\q)  in  the  V  norm,  uniformly  in  t  on  compact  subintervals  of  (0,oo). 
Indeed,  for  all  k  =  1, 2 . . . ,  we  obtain  §p;uN(t\ qN )  — >  Jpju(t;  q)  in  V. 

Since,  as  we  shall  see  in  the  next  section,  many  interesting  examples  involve  V  C  H2(il) 
with  Q  C  R1  or  f2  C  R2,  and  since  in  this  case  V  imbeds  continuously  in  (7(0),  the  above 
results  guarantee  pointwise  (in  t  and  x )  convergence  of  uN(qN )  and  all  its  time  deviatives 
to  u(q)  and  its  time  deviatives,  respectively.  Thus,  observation  operators  related  to  laser  vi- 
brometers  ( ut )  and  accelerometers  ( utt )  are  included  in  the  convergence  and  method  stability 
theory  for  output  least  squares  inverse  problems. 

Similar  results  are  available  for  the  nonhomogeneous  equation  (2.8)  (or  (2.10))  if  one 
places  appropriate  regularity  conditions  on  /.  To  obtain  these,  one  considers  the  representa¬ 
tion  theorems  (mild  solutions) 

wN(t,  qN )  =  TN(t ;  qN)PNw0  +  So  TN(t  -  *■  qN)PNF(s,  q)ds 

™(*;?)  =  T{t\q)w o  +  /o  T(<  -  S]q)F(s,q)ds 

and  uses  regularity  conditions  on  mild  solutions  in  H  =  V  x  H.  For  example,  if  f(-,q )  is 
in  C'1([0,  T],  H),  one  can  differentiate  once,  if  /(•,  q)  is  in  C'2([0,  T],  H ),  one  can  differentiate 
twice,  etc.  (see  [P]).  Indeed,  since  we  have  analytic  semigroups,  we  can  weaken  the  conditions 
on  /  needed  for  this  procedure  by  considering  special  regularity  theorems  for  mild  solutions 
of  Cauchy  initial  value  problems  in  the  case  of  analyticity  (e.g.,  see  Chapter  4.3  of  [P]).  For 
example,  it  suffices  to  have  f(-,q)  G  L1(0yT\H)  and  locally  Holder  continuous  in  order  to 
differentiate  wN  and  w  once  in  the  above  representations. 

4.  Examples 

In  this  section  we  present  briefly  several  examples  to  which  the  above  theory  can  be 
readily  applied.  In  some  cases  this  sharpens  the  convergence  and  stability  results  currently 
available  in  the  literature;  in  other  cases  it  provides  new  results  for  the  associated  inverse 
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problems. 


Example  4.1.  We  consider  a  cantilevered  Euler- Bernoulli  beam  with  Kelvin- Voigt  damping 
(i.e.,  stress  proportional  to  a  linear  combination  of  the  strain  and  the  strain  rate).  The  beam 
is  assumed  fixed  at  x  =  0  and  free  at  x  =  t,  with  the  transverse  displacement  at  time  t  and 
position  x  given  by  u(t,x).  Typical  observations  consist  of  acceleration  uH(t,Xj)  or  velocity 
ut(t,Xj)  at  several  locations  Xj.  Balance  of  forces  and  moments  yield  the  following  system 
(we  assume  the  linear  mass  density  p  is  normalized  to  unity): 


d\  a2  fE1Ta2u  T  a3u  ,  ,,  . 

at2  +  dx^EIdx2  +  CdI dx2di *  _ 


0  <  X  <  l, 


(4.1) 


r\ 

«(<»o)  =  ^(<,o)  = 0 


r  —  d2u  d3u 

+  CdI  a^8t 'I=<  “ 0> 


a3u 

dx2dt 


}U 


=  0. 


The  parameters  to  be  estimated  in  typical  examples  (see  [BWIC],  [BFWIC],  [Bln])  are  the 
stiffness  and  damping  coefficients,  El  and  cpl,  respectively.  For  a  parameter  set  we  choose 


Q  compact  in 


Qv  =  {q  =  (El,  col) :  q  <E  L°°(0,i)  x  L°°(Q,t),  EI(x)  >u>  0,cDI(x)  >  i ,}. 


For  the  state  spaces  H  and  V  we  take  H  =  H°(0,£),  V  =  Hl(0,l)  =  {<f>  €  H2(0,£)  :  <t>( 0)  = 
<£'(0)  =  0}.  Then  in  writing  the  system  (4.1)  in  the  form  (2.8)  we  take 

cx(q)(<l>A)=<EID2cj>,D V  >o 

<72(<l)(<t>A)  =  <  CDID24>,  D2\})  >0 

where  D 2  —  and  <  >0  is  the  inner  product  in  H°(0,£)  =  L2( 0,£).  For  real  valued  <j> 

in  V  we  have  at  once 

<Ti(q)(<t> >  4>)  >  v\D24>\o  >  cx\<j)\v 
which,  as  we  have  noted,  implies 


Re  a,(9)(0,  <f>)  >  rj^l2, 


for  complex  valued  <f>  in  V.  It  is  equally  trivial  to  argue  that  conditions  (A)  and  (C)  of 
Section  2  hold  for  this  o,i(q),o-a(q).  For  Galerkin  schemes  satisfying  condition  (Cl)  we  thus 
have  Theorem  3.2  applicable,  and  hence  parameter  convergence  and  method  stability  hold 
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for  the  least  squares  problems  {V)  and  (VN)  whenever  one  formulates  these  with  pointwise 
observations  of  either  velocity  or  acceleration. 


Example  4.2  For  this  example  we  return  to  the  so-called  “RPL  experiment”  discussed  in 
some  detail  in  [BGRW].  The  focus  of  our  attention  is  a  cantilevered  Euler- Bernoulli  beam 
with  a  flexible  gas  hose  and  thruster  nozzle  attached  to  the  free  end  as  depicted  in  Figure 
2.2  of  [BGRW].  The  structure  is  modeled  as  a  uniform  cantilevered  beam  with  Kelvin- Voigt 
internal  damping  and  tip  mass  with  a  mass-spring-dashpot  assembly  attached  at  the  tip. 
Along  with  the  usual  damped  Euler- Bernoulli  beam  equation 


.  .  d2u  d*u 

(4'2)  PW  +  EId7>  +  CDld7di 

we  have  the  force  balance  equation  at  the  tip 


=  0,  0  <  x  <  l, 


^  r  92u  T  54u 

(4.3)  [mr-^  -  cDIg^gt  -  El-^Ut 

ft'll 

=  CH{y(t )  -  +  kH(y(t)  -  u(t,i))  +  f(t) 

and  the  hose  assembly  state  equation 

dv, 

(4.4)  rnHy{t )  +  cH(y(t)  -  t , £))  -f  kH(y(t)  -  u(t,  £))  =  0. 

Here  p,  El,  ajid  cqI  are  the  usual  beam  parameters  whereas  mr  represents  the  tip  mass, 
mu  is  the  hose  mass,  kn  is  the  hose  stiffness  (the  “spring”  or  restoring  force  constant),  c# 
is  the  hose  damping  coefficient  and  f(t)  represents  an  externally  applied  force  at  the  tip 
(firing  of  the  tip  mounted  thrusters).  Boundary  conditions  for  the  coupled  state  equations 
(4.2)  and  (4.4)  include  the  tip  force  balance  equation  (4.3),  the  tip  moment  balance  equation 
(assuming  that  the  hose  assembly  has  negligible  rotational  inertia) 

(4.5)  1E^+Cd/K%)*='  =  0' 

and  the  zero  displacement,  zero  slope  conditions  at  the  fixed  end  x  =  0 

(4.6)  u(t,0)  =  |^(t,0)  =  0. 

The  structure  is  assumed  initially  at  rest  so  that  initial  conditions  are  given  by 


“(O,®)  =  -^-(M)  =  0 

2/(0)  =  3/(0)  =  0. 
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The  parameters  to  be  estimated  using  accelerometer  observations  (see  [BGRW])  include 
q  =  (ttit,  EI,ci}I,m}{,CH,kj{),  which  is  to  be  chosen  from  a  compact  subset  Q  C  R\. 

To  write  the  system  (4.2)-(4.6)  in  weak  or  variational  form,  we  use  the  state  variable 
u(t)  =  (y(t),u(i,^),u(2,  •))  in  the  state  space  H  =  R2  x  H°(Q, l)  with  inner  product 

<  (C>  v,  <t>)>  (A,  /*>  VO  >h=  C*  +  w+  <  4>A  >o  • 

For  the  space  V  we  choose  V  =  {{C,vA)  £  H  :  <f>  €  H2(Q,£),  <£(0)  =  D4>( 0)  =  0,t 7  =  $(£)} 
with  inner  product 

< «,  m,  «,(•*,  m,  «  >v=  « -  wm  -  m)+  <  d**,  d  v  >„ . 

The  stiffness  and  damping  sesquilinear  forms  are  given  for  <£  =  ((,  $(1),  <f>),  =  (\,tp(£),ip) 

in  V  by 

»l(«)(M)  =  M(  -  #M)(A  -  *(/))  +  El  <  >o 

-  «/))(A  -  +  cDI  <  C2V  >o  • 

We  also  need  an  operator  A4(q)  €  £(H)  given  by 

M{q){C,nA)  =  {rnH(,mTT},p4>). 


This  operator  can  be  extended  to  C(V*)  in  an  obvious  manner.  Then  the  system  (4.2)-(4.6) 
can  be  written  in  variational  form  for  the  state  u(t)  €  V  to  satisfy  for  all  <f>  £E  V 

4>)  +  <^2 (g)(*«(0»  4>)  +  ^i(?)(«(0. 4>)  =<  F(0. 4>  > 

(4.7) 

6(0)  =  «t(0)  =  0, 


where  F(t)  =  (0,/(t),0).  Since,  for  mH,mT,p  positive,  the  operator  M(q)  is  invertible, 
equation  (4.7)  is  obviously  equivalent  to  an  equation  of  the  form  (2.8).  Thus  the  theory  of 
Sections  2  and  3  is  applicable  if  and  ct2  satisfy  the  requisite  hypotheses.  If  Q  is  bounded 
below  in  R\ ,  it  is  readily  seen  that  both  o\  and  a3  satisfy  conditions  (A),  (B),  and  (C).  For 
example,  we  see  immediately  for  $  =  (£,  <j>)  £  V 


*»(«)(#,*) = **«  -  my + c^iWio 


>  <*«<  -  my + idvi?} = c \K- 

Similar  arguments  hold  for  <7i(q). 

Application  of  the  theory  in  Sections  2  and  3  substantially  sharpens  the  results  given  in 
[BGRW].  In  that  paper  the  main  results  (see  Lemma  3.1)  yield 

(4.8)  /  |u£(f;  qN )  -  Utt(t\ q)\Hdt  -*  0 

Jo 


£ 
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as  N  — >  oo.  The  arguments  are  rather  tedious  and  require  the  assumption  ii(q)  £  H2{ 0,  T\  V) 
on  the  limit  function.  Note  also  that  (4.8)  would  require  continuous  time  acceleration  ob¬ 
servations  be  used  in  the  least  squares  criterion. 

In  contrast,  the  theory  of  this  paper  yields  (since  H 2  embeds  compactly  in  C)  uN(t,x;  qN)  —* 
u(i,x;q),  uN(t,x]  qN)  —>  u(tyx\q)  uniformly  in  x  £  [0,^]  for  each  t  £  [0,  T1]  as  well  as 
iiu(t-,qN)  — >  utt(t]qN )  in  V  for  each  t,  which  permits  sampled  time  acceleration  observa¬ 
tions.  The  arguments  are  simple  (given  the  theory  developed  above)  and  do  not  require  the 
a  priori  smoothness  assumption  on  u (q). 


Example  4.3  As  a  final  example,  we  briefly  describe  the  models  for  two  dimensional  grid 
structures  (“plates  with  holes”)  developed  and  investigated  in  [R]  and  [BR] .  We  use  Love- 
Kirchoff  plate  theory  with  Kelvin- Voigt  damping.  The  “plate”  is  rectangular,  e.g.,  on  (x,y)  £ 
[0,£i]  X  [0,^])  with  rectangular  holes  periodically  placed  to  produce  a  thin  planar  grid.  It  is 
assumed  to  be  hanging  vertically,  clamped  at  the  top,  with  the  other  three  edges  free.  The 
basic  equation  for  transverse  displacements  w(t,x,y )  is  given  by 


d2w  d2Mx  nd2Mxy  82M *  , 

^  dt2  ^  dx2  dxdy  dy 2 

where  h  is  the  thickness  of  the  plate,  the  bending  moments  Mx,  Mv  about  the  x  and  y  axes, 
respectively,  are  given  by 

El  rd2w  d2w ,  cpl  ,  d3w  d3w 
l-v2Xdx2  dy2*  l-v2Xdx2dt  dy2dV 


Mv  = 


El  . d2w  d2w .  cpI  r  cPw  d^w 
I ~^2X'dy2  +  +  +  "d^di*’ 


and  the  twisting  moment  is  given  by 


Mxy  = 


El  d2w  d2w  cDI  d3w 

1  —  v2  dxdy  dxdy  1  —  v2  dxdydt 


d3w  . 
V dxdydt 


Here  v  is  Poisson’s  ratio.  For  the  boundary  conditions  at  the  top  (clamped)  we  have  the 
essential  boundary  conditions  w  =  =  0  (the  top  corresponds  to  the  x— axis).  The  plate 

is  free  on  the  other  outer  edges  and  on  the  edges  of  the  holes,  where  the  natural  boundary 
conditions  of  zero  moment  and  zero  shear  are  required.  For  example,  on  a  free  edge  parallel 
to  the  3/ — aocis  this  results  in  the  conditions 


Mx  =  0, 


dMx  dM *v 
dx  +  dy 


=  0. 


As  shown  in  [R],  [BR],  one  can  readily  define  the  corresponding  stiffness  and  damping 
sesquilinear  forms  <7i,<72  on  V  x  V)  where  V  =  {<f>  £  H2(Cl)  :  <f>  =  =  0  along  y  =  0},  Q 
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is  [0,/x]  x  [0,  -£a]  less  the  holes,  and  argue  the  needed  V—  coerciveness  along  with  conditions 
(A)  and  (C).  In  this  case,  the  state  space  is  H  ~  The  theory  of  Sections  3  and  4  can 

thus  be  shown  to  hold  for  this  example. 
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